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Introduction. The paper is devoted to the study on the three- 
dimensional model of transport and suspension sedimentation 
in the coastal area due to changes in the bottom relief. The 
model considers the following processes: advective transfer 
caused by the aquatic medium motion, micro-turbulent diffu- 
sion, and gravity sedimentation of suspended particles, as well 
as the bottom geometry variation caused by the particle set- 
tling or bottom sediment rising. The work objective was to 
conduct an analytical study of the correctness of the initial- 
boundary value problem corresponding to the constructed 
model. 

Materials and Methods. The change in the bottom relief aids 
in solution to the initial-boundary value problem for a parabol- 
ic equation with the lowest derivatives in a domain whose 
geometry depends on the desired function of the solution, 
which in general leads to a nonlinear formulation of the prob- 
lem. The model is linearized on the time grid due to the “‘freez- 
ing” of the bottom relief within a single step in time and the 
subsequent recalculation of the bottom surface function on the 
basis of the changed function of the suspension concentration, 
as well as a possible change in the velocity vector of the aquat- 
ic medium. 

Research Results. For the linearized problem, a quadratic 
functional is constructed, and the uniqueness of the solution to 
the corresponding initial boundary value problem is proved 
within the limits of an unspecified time step. On the basis of 
the quadratic functional transformation, we obtain a prior es- 
timate of the solution norm in the functional space L2 as a 
function of the integral time estimates of the right side, and the 


initial condition. Thus, the stability of the solution to the initial 


Beedenue. Hactosmuaa padota nocBaljeHa uccieqOBaHHIO 
IIPOCTpaHCTBeHHO-TpexXMepHOH MOJeIM TpaHcnopta Hu oca- 
%CHHA B3BECH B MIpHOpexHOM 30He Cc y4eTOM H3MeHeHHA 
pembeda gua. Moyen yuuTbIBaeT cilefyIoulwe MpoLeccsl: 
aBCKTHBHBIM WepeHoc, OOycNOBICHHLIM JBMKeHHeM BOJHOM 
cpebI, MHKpoTypOyeHTHYy!IO DHppy3u1o U rpaBuTaljHoHHoe 
ocaxeHHe 4aCTHI, B3BeCH, a TaKxKe H3MeCHeHHe reoMeTpHU 
Ha, BbI3BAHHOe OCaKJCHHeEM YaCTHI B3BeCCH WIM MObeEMOM 
YaCTHL JOHHBIX OTIOKeHHH. 

Lenvio pabombi AgiAioc’ TpoBeyeHHe aHayIMTH4ecKOrTO HC- 
Cyl¢OBAaHHA KOPpeKTHOCTH Ha4yasIbHO-KpaeBO 3aa4u, COOT- 
BeTCTBYHOMeH NOCTpoeHHOH MoJesM. 

Mamepuaazvi u memoovi. U3menenne pembeda Ha IpUBOHT K 
HeOOXOJMMOCTH pellaTb HavasIbHO-KpaeBylo 3ayayy WA 
ypaBHeHua Mapadommueckoro TuMa C MilaqWIMMH TpoH3BoL- 
HBIMH B OOJIaCTH, TEOMeTPHA KOTOPOM 3aBHCHT OT HCKOMOIi 
(PYHKUMM pelieHHsa, YTO MpHBOMT, B OOMeM cyIyyae, K HeJIH- 
HeliHOl NocTaHoBKe 3aya4uv. BamojHeHa JMHeapH3aljua MO- 
WelM Ha BPeCMeHHOHM CeTKe 3a CYeT «3AMOPaxKHBAHHA)) PelIbe- 
(pa WHa B Ipeyesax OAHOrO Wara 10 BpeMeHH MH Wociezyrouye- 
TO TMepecueta (yHKUMH MOBepxHOCTH Ha Ha OCHOBe H3Me- 
HMBUIelica (PYHKUMH KOHWeHTpal|HH B3BeLeCHHOTO BeLIeCTBa, 
a TaKoKe BO3MO2KHOrO H3MCHCHHA BeKTOpa CKOPOCTH JBHKe- 
HHA BOO cpesBl. 

Pe3ynomamei uccaedosanua. JA IMHeapv30BaHHOl 3aqauH 
MOCTpoeH KBaspaTHYHbIt PyHkWUMOHa WU 39HepreTHyecKHM 
MCTOJOM JOKa3aHa CJ[MHCTBCHHOCTb pelIeCHHA COOTBeETCTBYy- 
FoueH HauasIbHO-KpaeBOl 3aa4H B peyesax NPOW3BOJIbHOrO 
wlara io Bpemeuu. Ha ocHoBe npeoOpa30BaHHaA KBaypaTH4Ho- 
TO PyHKUMOHAaIAa NoWyyeHa allpHopHasd OL|eHKa HOpMBI pelte- 
HHA B (pyYHKIMOHAaIbHOM TpoctpanctBe L, B 3aBHCHMOCTH OT 
WHTerpasIbHBIX OWeHOK TO BpeMeHH TpaBoli yacTu, rpaHwy- 
HBIX YCJIOBHM HW HavasIbHOTO YyCJIOBHA, MH, TaKHM oOOpa3om, 
Woka3aHa YCTOMYMBOCTb pellleHHA HCXONHOM 3aa4H Ip v3- 
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problem from the change of the initial and boundary condi- 
tions, the right-hand side function, is established. 

Discussion and Conclusions. The model can be of value for 
predicting the spread of contaminants and changes in the bot- 
tom topography, both under an anthropogenic impact and due 
to the natural processes in the coastal area. 


Keywords: coastal systems, mathematical model, diffusion- 
convection problems of suspension sedimentation, bottom 
relief change, uniqueness of solution, and stability of initial- 


MeHeHHH HadasIbHOrO HM rpaHH4HbIx ycuoBui, pyHKUAN mpa- 
BoM 4acTH. 

O6cyarcoenue u 3axmouenua. Mojedb MoxkeT IpeCTaBsATb 
UCHHOCTh pH Mporno3se pacipoctpaHeHHA 3arps3HeHH u 
W3MeHeCHHA pelbeda AHa, KaK pH aHTpONOreHHOM Bo3zeii- 
CTBHM, Tak  B CHJIy ECTeECTBeHHO MIpoTeKaroulMx MpHpOAHBIX 


Hpoweccos B npHOpexkHoli 30He. 


Kunrouespble cs10Ba: TIpHOpexKHBIe CHCTCMBbI, MaTCMaTH4eCKaA 
MOJI€JIb, 3a,}adn Tuddy3uu-KOHBeKLMH OCMKTCHHA B3BCLICH- 
HOrO BeIICCTBa, H3MCHCHHe pembeda 7jHa, CXHHCTBCHHOCTB 


boundary value problem. pelleHHa HW yCTOM4MBOCTL HavaJIbHO-KpaeBOH 3aa4H. 
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Introduction. The aquatic habitat protection [1—2] is one of the most important factors that determine the inte- 
grated research development of the coastal areas. Damage control over the natural processes, such as pollution, sedi- 
mentation, and depletion of water areas, leads to necessity for studying all aspects that affect changes in coastal waters. 
Maintenance of water bodies in proper condition and timely intervention in its operation mode is directly related to the 
increase in port capacity and the efficient development of the coastal infrastructure (ensuring an accessway to the berths 
of ships with a low landing; desilting and aquatic vegetation clearing of the coastal strip; etc.) [3-5]. As a rule, research 
practice in this field requires the construction of mathematical models that are as close as possible to real processes 
[6-11]. 

A continuous mathematical model describing spatial-three-dimensional processes associated with transport and 
gravitational suspension sedimentation in the aquatic medium with varying bottom relief is considered. This model 
takes into account micro-turbulent diffusion and advective transfer of suspensions, the effect of gravity on suspension, 
the presence of the bottom and a free surface, and a bottom contour variation. 

The suspension transport model enables to study the hydrophysical processes of aquatic systems, to predict the 
dynamics of the bottom surface change based on the description of the lifting, transport, sedimentation, changes in the 
concentration of suspension [12—13]. The uniqueness of the solution to the corresponding initial-boundary value prob- 
lem is proved, and a prior estimate of the solution norm is obtained depending on the integral estimation of the right- 
hand side, boundary conditions, and the initial condition. 

Materials and Methods. Continuous 3D model of suspension diffusion-convection and the corresponding initial 
boundary value problem. Consider a continuous mathematical model of sediment spreading in the aqueous media 
considering diffusion and convection of suspension, gravity action on suspension, presence of the bottom and a free 
surface. We will use Oxyz Cartesian coordinate system where Ox axis passes along the nonperturbed water surface 


and is directed toward the sea, and Oz axis is directed vertically downwards. Assume that A=H+n is the total water 


depth, m; H is depth with undisturbed water surface, m; n is elevation of the free surface relative to the geoid (sea lev- 
el), m (Fig. 1). 








— ee - gee 


fH y) 





Fig. 1. Introduction of Oxyz coordinate system Oxyz 
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Suppose that in G = {0 Sesh Us ysl Uses 8 (x, y) closure region, there are suspensions which have 





c=c( x y,z,t) concentration at xe y,z) point and at ¢ time, mg/l; ¢ is temporary variable, sec. We will also use 
L.= max  H(x,y) notation. 


Osx<L,, 0SVSL, 


The behavior of the suspended particles will be described by the following system of equations: 


Ge , (ue) (ve) Olle. )e) 7 (2 F ee @ (» 20) 





ot =x oy Oz ao ay) a” & (1) 
oH € 

—=--—w,c, 

ot p * 


where wu, v, w are components of U fluid velocity, m/s; wg, is hydraulic size or sedimentation rate, m/s; UW, 1, are coef- 
ficients of the horizontal and vertical turbulent diffusion of particles, respectively, m’/s; F is power of particle sources; ¢ 
is porosity of bottom materials. 

Summands on the left side (except for the time derivative) of the first equation of the system (1) describe the 
advective particle transport due to the inertial motion of the aqueous media, as well as sedimentation under the action of 
gravity. The summands on the right side describe the suspension diffusion. The vertical diffusion coefficient is chosen 
different from the horizontal diffusion coefficient due to the fact that the effect of difference between these coefficients 
is often observed in various media and can be caused by various factors. 

As G region, we consider ABCDAOC,_D, “parallelepiped” “skewed” to the shore, whose 4,OC,D, upper base 


lies on (z=0) free surface, and (z = H(x,y)) part of the bottom surface is its lower base. Suppose S is G surface, 7i 


is the outward normal to the surface of the “skewed parallelepiped”. We assume the given U" as the fluid velocity on 
G side surfaces. Complete with the boundary conditions of first kind for the particle concentration function, this allows 
determining the suspension flow both towards the coast and along the coast (Fig. 2). 


Cc, Coast 
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Fig. 2. Solution area for suspension transport 


Add the initial and boundary conditions (assuming that the sedimentation is irreversible) to the system (1) 
As the initial conditions at ¢ = 0 time, we accept 


c(x, y, 2,0) = cy (x,¥,2)5 (2) 
H(x,y,0) =H, (x,y). (3) 
We set boundary conditions on ABCDA,OC,D, faces (we set suspended flows both towards the coast and 


along the coast): 
- on the faces S,=AAOB (x=0, OS y<L,, 0<z<L,) , S,=AADD (y=L,, 0S xSL,, 0Sz<L,) and 





S,=BOCC (y=0, 0<x<L,, 0Sz<L,) 








e=e¢ ,me¢ =c (x,9,2,t), fe[0,7]; (4) 
- on the faces S, =DD.C,\C (¥=L, 0<y<L, 0<z<L,) and S, = AOC D, (2=0, 02381, O<y<L) 
c=0; (5) 


- on the surface S, = ABCD (2=H(x,y,t), O<x<L, O<y<L,) 
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Os FA eg Oe ae, (6) 
on Lu, Oz L, 


The boundary condition (5) occurs with a relatively small slope of the bottom: 


(4) dH \ 
max ,/} —-| +} —]|] <1. 
Ss Ox Oy 


The following condition of the solution domain nondegeneracy is set up for all (x, y,t) at which the initial 


boundary value problem is formulated: 

H(x,y,t)>hy =const>0, O<t<T. (7) 
When studying combined models of sediment and suspension transport, it is possible to increase the concentra- 
tion of suspended particles in the bottom layer due to the rising bottom sediment particles if the shear stress exceeds of a 
certain critical value is exceeded[13—16]. Then, instead of the boundary condition (6), we will consider the boundary 

condition of the following form 
Oc | 
7 
Linearization of the initial-boundary value problem of transport and suspension sedimentation. To create a line- 
arized model on 0<t<T time interval, we construct a uniform grid wt with a step t, that is, a set of points 


oe =11, = 01, n= 01,4, NST, 


ac, a&=const > 0. (8) 








Fig. 3.Construction of time grid 


ce (x, y,z,t,,)and H” (x, y,t,,) functions are determined at each step of , time grid. If n=1, then func- 
tions of the initial condition will suffice to c"? (2,9, 25t, )s H” (x,y.t,) , viz ol) (x,y,z,0)=c,(x,y,z) , 
H" (x, y,t,) = H, (x,y) respectively. But ifn =2,...,N , then ec” (x,y,z,t,,)=cl" (x, y,z,4,,) functions are assumed 
to be known, since the problem (1)—(6) for the previous t,_, <¢<t,, time interval is supposed to be solved. 


We write the system (1) on t,, <¢<t, interval in the form: 


ac) A(uc) ave) a((w+aw, Je) (a — o( sc) ; 
=H, + + 36 +4, 





~ + ~ 


Ot Ox oy Oz ax? oy’ Oz Oz (9) 
oH") € (n) 
—— =-—w,c 
ot p ° 
and complete it with the initial conditions: 
ol) Coed, =e, (x,y,z), (*,9,2,4,4)= ol) (2,9, 2,4 )s #=2,4,N. (10) 
H") (x, vety) =H, (x,y), (ort) =H") (54 ), n=2,...,N. (11) 


The boundary conditions (4)—(6) are assumed to be fulfilled for all ¢,, <¢<¢, time intervals. 


oe 
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By defining c'”(x,y,z,t,,)=cl"(x,y,z,t,,) function on t,,<¢<t, time interval, we can find 


n-l 
H ae y; t,) function. For this end, we integrate both members of the second equation of the system (9) over 


t_, <t<t, variable. We will get 


oo 





6 oH) 
| ~ dt=—~w, | cat. (12) 
p 


fa 


From the equality (12), it is not difficult to get 


n a~ é ul ‘ n 
H”) =H! 7 Mey [era (13) 


n=l Heat 


ml 


We introduce G,_,= {0 4241 Nagel 0e2eH™ (ait ) domain at each ¢,, <¢<t, time step. 


We have a chain of linear initial-boundary value problems for each time layer, where the system of the type 





+ ~ + —+—- |+—] 
Ot Ox oy Oz ox oy Oz Oz 


(xy.z)EeG_., G 


n-1 


ey fl 
= dt, n=1,2,...,N. 
ovate (15) 


ac”) Aue”) a(ve) a((w+w, je) [ae a) o( sc) ; (14) 
— Hi, y +h, 


{O<x<L,, 0<y<L,, 0<z<H""(x,y,t,,)}, 


HH =H") 


is considered for ¢,, <¢<¢, interval with the initial conditions: 
ol”) (x, Vy Rina) =") (3% y,Z,t,_, ), (16) 
H” (4.9:2,4) =H") (ithe ) : (17) 
Note that at each time step, the boundary surfaces will change (except S, face). Considering ¢,, <¢<¢t, time 
interval, we set the boundary conditions on the edges of G,_, domain: 


n-l 


- on S 





1jn-l 


(P=0, 0S yeh, OSesn"(0y2,)) g Sia(pols O<eel, 0x2eR (2 i.4,,)) “and 
S,,.(y=0, OSxSL,, 0Sz<H"(x,0,t,,)) faces 
oc” =e" ,rme c’ =c"(x,y,z,t), tet, .t,]; (18) 
-on §,,,(x=L,, OS ySL,, OSz<H"(L,,y,t,,)) w S,,,(2=0, OSX<SL,, OS ySL,)=AOCD, faces 
ees (19) 
-on S,,,(2=H"" (x,y,t,,), OS XSL, 0S ySL,) surface 











(n) (n) i“ 
a rn oes (20) 
on LL, oz L, 
The boundary condition (8) will be replaced by the following 
(») 
cae ac"), a =const > 0. (21) 
Oz 

Thus, it is supposed that the bottom relief within this time step, when calculating the distribution of suspension 
concentrations, does not change and is taken from the previous time layer. First of all, at this ¢,, <<, time step, the 








initial-boundary value problem for the convection-diffusion equation (14) with H'” fixed bottom relief function is 


solved, and only then the update (recomputation) of H'”’ relief function is performed in accordance with the equality 
(15). 

The determination of the conditions of existence, uniqueness and continuous dependence of the solution on the 
input problem data is carried out on a fixed time layer under these assumptions and subject to the condition (7). 

The authors do not plan to study the existence of solutions to the initial-boundary value problems (14) - (20) 
and (14) - (19), (21) in this paper. Questions of the existence of solutions to the initial-boundary value problems for 
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parabolic equations with lower derivatives (diffusion-convection equations) are considered, for example, in the mono- 
graphs [17-18]. 

Research Results. Investigating uniqueness of the solution to the initial-boundary problem of suspension 
transport. 
Consider the initial boundary value problem (14) - (20) formulated for the arbitrary t, ,<¢<t, time layer. 


Multiply the left and right member of equation (14) by c” function and get: 


























)  A(uce™)  A(ve™) a((w+w, )c 2() a2 (n) (») 
alan ( ), ( ), (( :) ) =p,c” eZ a cae LL, Ose R, (22) 
ot ox oy OZ ox oy” OZ Oz 
The left member of the equality (22) can be transformed as follows: 
2 
(n) A(uc”) dl(ve) a((w+w, Je” a(c” _ 
ol”) oc +c”) ( le ( ) (( | _ 1 ( ediv(c7) = 
Ot Ox oy Oz 2 Ot 
(23) 
a(cl”) 1 Bon 
ee ( +5aiv((c) 0) , 
2 oat 2 
where U =|lu, v, w+ w, [’. 
With regard to (23), the equation (22) will be written as 
2 
a(c” - ral) a2_(n) (n) 
. ( + 5aiv((e) 0 = p,c” = + g = +o) eC Ly & +0 F, (24) 
2 oat 2 ox oy Oz Oz 


Then we integrate both members of the equation (24) over ¢,, <¢<¢, interval, and, after that, over the spatial 


variables in G 


, domain. In the first term, the order of integration is changed due to the Fubini theorem [19]. We obtain 


i (n) : t, 
We ae G,.+] 4 aera. 
Gua 2 tes at ta 2 Gra 


(iter (222-2 Joc,. Jue Her Sn 22 laoa ee 


{| (ere, 


ta\ Gia 








The first term on the left side of the equation (25) is obviously equal to 


3] fA ade |(6F or.24)-F toreted) eee re 


Next, we turn to the transformation of the second term of the left-hand side of the equality (25). Considering 
the Gauss-Ostrogradsky formula and the boundary conditions (18) - (20), it can be written as [20]: 


j ¢ J | J div((e” U)aG,, yn = | Ci (c’) (v, av + a i J | ey dt + 
+f eV (e se sig (cy (Oviatt - -si(u (cy nyt - on 


byt S35 1 


oe I [Hey vdsde o£ iL J J | (c’) vied es If J | (ce) want dt. 


28 s\ Ses 


where U" is the known velocity of the aquatic medium on the faces where the boundary conditions of the first kind 


are specified; in fact, these are all side faces, except for S. and S$ 


sh top cover on which the suspension concentra- 


5,n-1 
tion is zero, and therefore the flows through them are zero. 
Let us turn to the transformation of the right side of the equation (25). The following equality occurs 
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fsa alee) 


() (») (x) 
! Ox Ox oy oy Oz Oz 


- ac") 2 ac") 7 ae”) 2 
- ——| + —— | +n, dG... 
I) LH, ax LH, dy H, a n-l 




















= (») ac” (») 6c") (n) ac” 5 < 
Suppose 0={0,,0,,0.} =4U,C yU,c sc, . Then, in virtue of the Gauss-Ostrogradsky 
: ° ox oy Oz 


theorem, we have: 








ac” al moc™ a(m a (n _ 
[f} om Ser Jou, 2 [eS] Sf cry, Ha, fad - 


ee dxdz + J O,.dydz + Il O,dxdz + I O.dydz + J QO,dxdy + ‘I QO, dxdy = (29) 
“I O, dxdz + ii O,dxdz + I O.dydz + I Ci 


Son 


Transforming each term of the dvbehand side of 09) subject to e conditions on the boundary (18) - (20), 


we obtain 
(n) (n) (n) 
i (a) {an Oc se - fe" “ }. - ( (n) Ll, “ } dG, = 
ox y y IZ IZ (30) 


cu die + cu Sy tet cu dda - Oh dxdy. 
= [fens Jen Mem Mvle) 


In virtue of 06), (28), i and (30), the cae (25) ie on form 


zllfle) re oat ey ety, & = 
‘als (c) v+c'n, agmleloe (ey | v ong me 
3)( geet ae of) (22) nfo. fe 
=> fffle I} (xy. 25t,) 14, f| fe, 


ta \ Gra 











(31) 











The identity (31) will be fundamental under studying the uniqueness and obtaining a prior estimate of the solu- 
tion norm of the initial boundary value problem (14) - (20). In case of replacing the boundary condition (20) with the 
boundary condition (21), the quadratic functional (31) changes as follows: 


allile "y( x, y,2,t, dG, jake c ‘ween, = 
(alyer v+c'u, oie pl[ ge *y—c', oat 
TU G-on le yao l() si te 


= fll?) (sa.2 ‘oot. ff Fag 


Suppose that the equation (14) with the same conditions (16) - (20) satisfy two different solutions to 





(32) 








C=C, (x, y; a); C, =C, (x y,z,t) problem. For their c= c, —¢, difference, the following initial-boundary problem is 


valid: 
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For c function, the equality (33) will take the form considering the equalities (34)—(36) 
fle X,y,Z,t, )dG, , += 2 [ff wan dt + 
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Since w, >0 and other known values under the sign of integrals are positive uh> 0, v> 0, then the equality 
(36) is satisfied only under the condition 
é(x, y,z,t) =0, (x,y.z)EG,, tt <tSt,, (38) 


which completes the proof of the uniqueness of the initial-boundary value problem (14) - (20) solution. 
In case of replacing the boundary condition (20) by the relation (21), instead of the expression (37), we obtain 
the following equality 
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We require the fulfillment of the inequality 


1 
oi aha 20, (GHzZES eas t  <tSt, 


or 


We 
as : (x,9,2) € S41 ti <t&t 
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then all the terms in the equation (39) are nonnegative, and zero equality is possible if and only if 





(40) 


n? 


y 


Ca; V5Z, t) =0, (x: ys z) eG,,, t,,<t<t,, that means the solution uniqueness and in this case. 

Reasoning is similarly repeated for all layers of w, time grid. The modification of the boundary conditions as- 
sociated with the continuous change in the bottom relief depending on the time variable requires additional study and is 
going beyond the scope of this article. 

Theorem. Suppose we are given a system of equations 
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(92)€0.., Oy af0exeb,0<7<b062 <1" (5) 
H® =H =n = f cd t, n=1,2,...,N 


elf, 


in 0, =G_x(6y<t<4), Ga= (0 <a 0<y<L0<2<H? "(x yt, ; )). simply connected domain with a 





sufficiently smooth boundary defined by the smoothness of z= H'"" (x,y), O<x<L,, 0<y<L, function with the 
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initial and boundary conditions (16) - (20). Let the functions of c")(x,y,z,t,,) solution, the velocity vector of 


Tr é a ees de cela aye F ’ 
u,v,w+w,|| aquatic medium, c"" (x, y,z,t,,) initial condition, F (x, y,z,t), right member of c (x,,z,t) bounda- 














ry condition, pb, =H, (as (x, y,z) eG,, coefficient of u, =H, (z), (x, y; z) € G,_, vertical turbulent exchange satisfy the 
following smoothness conditions: 
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<ec((o<x<1,, a ae ee z=H" (x,y))xf[t,, <¢<¢,]) , as well as c’(x,y,z,0)=c,(x,¥,z), 
(x,y,z) €S,,\(0<x<L,,0<y<L,, 2=H™ (xy); See, 


(0 <x<Lj0<y<L, z=H (9) (x, y)). conditions of consistency of the boundary and initial conditions, then the 


solution to this problem exists and is unique. 
Comment. In case of replacing the boundary condition (20) with the boundary condition (21), the inequality 
(40) should be added as a sufficient condition for the previous theorem. 

Studying the continuous dependence of the solutions to the initial-boundary value problem of suspension 
transport on the initial, boundary conditions and the right-hand side function. The next stage is connected with 
the study of the continuous solution dependence on the functions of the right-hand side, boundary and initial conditions 
for the system (14)-(15). 

Suppose that 
c’ >c, =const > 0, 
OSnk1, Use, O<2eH (4085) ta hte. @)) 
For convenience, we introduce the notations: union of all parts of the lateral cylindrical surface (boundaries of 
and the lower base of region— as G_,— S,_, . In virtue of the smoothness conditions 


n-l 


G,,.; region) is denoted as S 


en-1? 


listed under the above theorem, extrema of functions on the bounded closed sets are reached: 
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M,,., = max {Hy}, Myni =. max {lul.|y]}. Ms. = min {hn}. 


We will focus on the equation (31) if the boundary condition (20) is used, and on the equality (32) in case of 
the boundary condition (21). Evoking Friedrichs inequality, we have a chain of inequalities: 
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We turn to the equation (26) from which, in virtue of (42) and (43), we obtain the inequality: 
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The inequalities obtained imply the continuous dependence (stability) of the solution to the problem (14) - (20) 
on the functions of the initial condition, the boundary conditions and the right-hand side, in L, norm for any instant of 


0<T <+co time, and also in L, time-integral norm. 


Obviously, if the inequality (45) and the theorem condition are satisfied, the initial-boundary problem (14) - 
(19), (20) will also have a solution that depends continuously on the functions of the initial condition, the boundary 
conditions and the right-hand side in the corresponding norms. 

Discussion and Conclusions. Novelty of the proposed non-stationary spatial-three-dimensional mathematical 
model of suspension transport lies in the fact that, alongside with considering the processes of advective transfer, micro- 
turbulent diffusion and gravity sedimentation of suspended particles, the model describes the change in bottom geome- 
try caused by the particle settling or bottom sediment rising. 

The linearization of the corresponding initial-boundary problem on the time grid is carried out, and the conditions 
for the uniqueness of the solution to the initial-boundary problem and continuous dependence on the input data — on the 
functions of the initial condition, boundary conditions, and the right-hand side in L, Hilbert space norm in L, time 


integral norm for two variants of boundary conditions are obtained for the arbitrary ¢,, <¢<¢, time step. 
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